Refined gluino and squark pole masses beyond leading order 



O 
O 

(N 

wo: 
< 

(N 



> 

(N 
O 
oo 
O 

o 

Oh: 
a: 



X 



Stephen P. Martin 
Physics Department, Northern Illinois University, DeKalb IL 60115 USA 
and Fermi National Accelerator Laboratory, PO Box 500, Batavia IL 60510 

The physical pole and running masses of squarks and gluinos have recently been related at two-loop 
order in a mass-independent renormalization scheme. I propose a general method for improvement 
of such formulas, and argue that better accuracy results. The improved version gives an imaginary 
part of the pole mass that agrees exactly with the direct calculation of the physical width at next-to- 
leading order. I also find the leading three-loop contributions to the gluino pole mass in the case that 
squarks are heavier, using effective field theory and renormalization group methods. The efficacy 
of these improvements for the gluino and squarks is illustrated with numerical examples. Some 
necessary three-loop results for gauge coupling and fermion mass beta functions and pole masses in 
theories with more than one type of fermion representation, which are not directly accessible from 
the published literature, are presented in an Appendix. 
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I. Introduction 

The small ratio of the electroweak symmetry break- 
ing scale to the Planck mass can be stabilized [1] in 
softly-broken supersymmetric extensions of the Standard 
Model. This implies that all of the Standard Model par- 
ticles will have superpartners, which should be within 
reach of the pp Fermilab Tevatron collider or the pp Large 
Hadron Collider during the next few years. Most of the 
new parameters appearing in the Minimal Supersymmet- 
ric Standard Model (MSSM) [2] are the masses of the new 
superpartners and other supersymmetry-breaking cou- 
plings of positive mass dimension. Therefore, a detailed 
understanding of the MSSM Lagrangian is nearly synony- 
mous with an understanding of supersymmetry breaking. 

The fact that experimental observations of flavor viola- 
tion and CP violation are not in significant disagreement 
with the predictions of the Standard Model can be taken 
as indirect evidence for the existence of some powerful 
organizing principle governing supersymmetry breaking 
and its mediation to the MSSM sector. An especially in- 
teresting possibility is that the organizing principle can 



be discerned by running the parameters of the theory up 
to high energy scales using the renormalization group. To 
carry out this analysis, it will be crucial to relate phys- 
ically measured observables, especially the superpartner 
masses, to running parameters in the full theory defined 
by the non-decoupled Lagrangian that includes all of the 
superpartners. 

However, running masses are not the most direct ob- 
servables expected from collider experiments. In general, 
the mass defined by the position of the complex pole 
in the propagator is a gauge-invariant and renormaliza- 
tion scale-invariant quantity [3]- [11]. The pole mass does 
suffer from ambiguities [12] due to infrared physics as- 
sociated with the QCD confinement scale, but these are 
probably not large enough to cause a practical problem 
for strongly-interacting superpartners. The complex pole 
mass should be closely related in a calculable way to the 
kinematic observable mass and width reported by exper- 
iments [13]. 

It is often useful to calculate in on-shell schemes, in 
which some physical masses and other observables are 
used as input data and others are outputs. However, for 
the key purpose of unraveling the organizing principle 
behind the supersymmetry-breaking Lagrangian, this is 
not as directly useful. The MS scheme [14] can also be 
used, but it violates supersymmetry explicitly. Instead, 
it is preferable to use the DR scheme [15] (or the revised 
DR scheme [16], which removes the unphysical effects 
of epsilon-scalar masses in softly-broken supersymmetric 
models), with all superpartners non-decoupled. While it 
is difficult to know in detail what limitations on this pro- 
gram will follow from future experimental uncertainties, 
it seems clear that multi-loop calculations will be neces- 
sary to make the theoretical sources of error negligible. 

The one-loop relations between the superpartner pole 
masses and the running parameters in the MSSM La- 
grangian have been known for some time [17]- [19] . The 
calculation of the Higgs scalar boson masses in the MSSM 
has now advanced to include the important two-loop cor- 
rections (for some reviews of recent progress, see [20]- 
[24]), and even some three-loop corrections [20], using a 
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variety of different methods. Recent calculations have 
provided the supersymmctric QCD (SUSYQCD) two- 
loop corrections to the squark [25] and gluino [26, 27] 
masses. The quark masses in the MSSM arc known at 
two-loop order [28], [27]. More generally, refs. [25] and 
[27] provide the self-energy functions and pole masses 
for scalars and fermions, respectively, calculated in mass- 
independent renormalization schemes at two-loop order 
in any renormalizable field theory, in the approximation 
that vector bosons are treated as massless in the two- 
loop parts. This approximation is likely to be quite 
good for most applications to the MSSM, because the 
largest two-loop effects involving vector bosons come 
from SUSYQCD, and because the W and Z bosons are 
evidently lighter than most of the supcrpartners. 

It is important to consider the validity of the perturba- 
tive expansion in these results. Especially for the lightest 
Higgs boson, the squarks, and the gluino, the corrections 
that give the pole masses from the running masses turn 
out to be quite significant. As a prominent example, even 
the pure two-loop correction to the gluino mass (com- 
pared to the running mass evaluated at a renormalization 
group scale equal to itself) is of order 1-2% in the case 
that squarks and gluinos arc comparable in mass, and 
grows to about 5% for squark masses that are of order 5 
times heavier than the gluino. One would like some assur- 
ance that perturbation theory is really converging, and 
an estimate of the theoretical error. Unfortunately, the 
renormalization-scale dependence of these results is not 
a reliable error estimate; in particular, the scale depen- 
dence of one-loop corrections is routinely much smaller 
than the two-loop corrections when the latter are known. 

Another reason to be wary is the fact that calculations 
in mass-independent renormalization schemes like MS or 
DR use propagators with masses that can differ signifi- 
cantly from the physical ones. In many cases this is true 
for any reasonable choice of the renormalization scale Q. 
A troubling aspect of this is that the imaginary part of 
the complex pole squared mass, 
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M - iTM 



(1.1) 



can give a numerical value for the width T that differs 
quite badly from the physical width. It is not hard to 
find examples for which the tree-level masses are suffi- 
ciently different from the physical masses that a partic- 
ular contribution to T as computed from the complex 
pole mass is exactly (because the decay would be kine- 
matically forbidden if the particles had masses equal to 
the tree-level Lagrangian masses appearing in the prop- 
agators of the self-energy loop diagrams), while the true 
decay width contribution (computed directly from dia- 
grams with multi-particle final states, using an on-shcll 
scheme) is non-zero. Or, the reverse can happen. (I 
will show an example of each type in Figure 2 of sec- 
tion IV.) While the complex pole mass is in principle a 
gauge-invariant and renormalization-scale invariant ob- 
servable, this calls into question how well one can trust 
the perturbation theory that yields it in practice. 



These issues are general. However, in the MSSM, they 
are particularly acute for the squarks and the gluino, be- 
cause of their strong coupling. Furthermore, the LHC 
will quite likely produce gluinos and squarks in abun- 
dance if supersymmetry is correct. Therefore, I will 
use the squark and gluino SUSYQCD system within the 
MSSM as an example in this paper to show how to ame- 
liorate the problems mentioned above. First, in sec- 
tion 2, I discuss how to reorganize the results of per- 
turbation theory by expanding tree-level masses around 
physical masses in the loop corrections obtained in mass- 
independent (MS or DR ) schemes. In section 3, 1 present 
a result for the three-loop corrections to the gluino mass, 
valid in the limit that squarks are treated as nearly de- 
generate and much heavier than the gluino. This is the 
case where one might expect three-loop and even higher- 
order corrections to be most dangerous, but I show that 
they are actually under good control, and can be tamed 
by using effective field theory and renormalization group 
methods. Section 4 displays some numerical results show- 
ing the efficacy of these improvements. In an Appendix, 
I present some necessary results for three-loop contribu- 
tions to fermion mass beta functions and pole masses in 
(non-supersymmetric) theories with fermions in distinct 
representations. 



II. Improved two-loop pole mass results 

In general, the two-loop order expression for the pole 
mass of a particle can be computed from knowledge of 
the self-energy function, 



1 n (2)fc + 

(16tt 2 ) 2 j 



(2.1) 



Here s is the external momentum invariant, the super- 
script in parentheses indicates the loop order, and the 
indices j, k indicate different particles with the same 
quantum numbers, which in general can mix. For 
fermions, (s) can be assembled from separate chirality- 
preserving and chirality-violating self-energy functions, 
as described in section II. C of ref. [27]. Then the 
gauge-invariant and renormalization-scale invariant pole 
squared masses can be defined formally as the solutions 
to the equation 



Det[(s-m 2 )^-n^(s)] =0, 



(2.2) 



where mf are the tree-level diagonalized squared masses. 
However, because Hj should be interpreted as a complex- 
valued function of a real variable s, this equation must 
be solved by first expanding the self-energy function as 
a series about a point on the real s-axis. In evaluating 
the loop integrals in the self-energy, s is given a posi- 
tive infinitesimal imaginary part, while the complex pole 
squared mass solution [see eq. (1.1)] always has a non- 
positive imaginary part. A related subtlety is that when 
the particle in question has couplings to massless gauge 
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bosons, terms of a given loop order in the self-energy have 
branch-cut singularities (except when the Fried- Ycnnic 
gauge-fixing condition is used [10, 25]). 

The most straightforward way to obtain the pole mass 



at two-loop order in a mass-independent renormalization 
scheme is to first expand 11^ (s) in a series about the tree- 
level squared masses. Define, for a generic squared mass 
to 2 : 



nJ 1)fe (TO 2 ) ee hm n5 1)fe (s), (2.3) 

(2.4) 



nf b '(TO 2 ) EE lim 



where the self-energy functions on the right-hand side are computed in a mass-independent renormalization scheme, 
and no sum on j is implied in cq. (2.4). Then, working consistently to two-loop order, the pole mass for the particle 
with tree-level squared mass to 2 is 

M]-iT jMj = ^ + ^n?V?) + 7^{n^^ (2.5) 

obtained by a perturbative' 1 ' expansion of eq. (1.1). 

The previous expression is gauge- invariant, and renormalization scale invariant up to terms of three- loop order. In 
its application to the squark and gluino masses in the MSSM, this approach has the advantage of depending only on 
tree-level running parameters, so that iteration is not necessary if they are taken as given. However, as remarked in 
the Introduction, the use of tree-level running masses in propagators is problematic at least for the imaginary part 
of the pole mass, which arises from the absorptive part of the self-energy functions. If the tree-level masses differ 
significantly from the physical masses of the particles, then the kinematics of the self-energy functions will poorly 
reflect the actual kinematics giving rise to the physical width of the particle. This can lead to a non-zero width 
when there should be none, or vice versa. In general, one may care more about the real part of the pole mass, but 
intuitively one cannot expect the real part to be very accurate if the kinematics in the loop integrations poorly reflects 
the physical particle masses, and if the imaginary part is completely wrong. 

To improve the situation, let us reorganize the previous result by expanding all tree-level squared masses appearing 
in the II functions in a series about the real parts of their respective pole squared masses. (Note that although I 
have written the II functions as depending on a single external squared mass argument that replaced s, there is also 
dependence on the internal propagator masses which is not explicitly indicated.) Doing this, one arrives at: 

M]-ir jMj = m| + ^n^(M 2 ) + ^{nfV^ 

-^M^ ) \Ml)] 1 ^Tif\M])}, (2.6) 

The second sum over k is taken over all fermions and bosons in the theory that couple to particle j, not just those 
that can mix with j. The II in eq. (2.6) are defined to have the same functional dependence on the real parts of 
the pole squared masses (both external and internal) as the functions II in cq. (2.5) did on the tree-level squared 
masses. Because the one-loop and two-loop parts of eq. (2.5) are each separately gauge-fixing invariant, it is clear that 
eq. (2.6) is also independent of gauge-fixing. Formally, eqs. (2.5) and (2.6) are equivalent up to terms of three-loop 
order. However, in eq. (2.6), all kinematic dependences of loop integrals on the right-hand side correspond to the 
physical masses (the real parts of the pole masses). I therefore expect that, when it makes much of a difference, 
eq. (2.6) should be more accurate than eq. (2.5). If one starts with the Lagrangian parameters as input, evaluation 



t Here the loop-induced mixing between particles j and k is as- 
sumed small compared to the tree-level squared mass splitting. 
Otherwise, one must perform almost-degenerate perturbation 
theory, by expanding around a modified tree-level Lagrangian 
designed to minimize the one-loop mixing in that sector. This 
could plausibly occur for Higgsino-likc ncutralinos in the MSSM. 
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of the pole masses will require an iterative procedure involving all of the particle masses simultaneously, which in 
a general case at two-loop order could take a significant computation time. On the other hand, if physical masses 
are taken as inputs, eq. (2.6) still requires knowledge of the tree-level couplings and mixing matrices of the theory. 
The one-loop functions are always known analytically in terms of logarithms, so taking derivatives of them poses no 
technical difficulties. The two-loop functions often cannot be computed analytically, but computer codes such as TSIL 
[29] provide for their numerical computation. 

To illustrate the method, I choose here to consider the two-loop corrections to the gluino and squark pole masses, 
for simplicity including only SUSYQCD corrections and ignoring the effects of squark mixing, quark masses, and 
electroweak effects. For this system, the realization of eq. (2.5) is: 



Ml - M/s. = m| . + hCgU^ (m| . , m?) + h 2 C q 



C g ng-)(m| j ,m? ) + ^nf)(m| j ,m? ) 



Qj' 1 



(2.7) 



M?-iT- g M g = m l + h[c A nf a \ m l) + I q Y,^ b \rnl im l^ 



+CaI q Y: b) (ml , ml ) + C q I q J2 fi? c > (m? , m J ) + I 2 q £ £ fif d > (m? , m% , m% k )] . (2. 



Here, to? and m 2 ~ are the tree- level DR squared masses. In many references, m g is written as M3, but in the present 

paper, capital letters are reserved for pole masses and lowercase letters for running masses. The strong gauge coupling 
appears in the combination: 



h = . 9 |/16tt 2 = a s /An. 



(2.9) 



The indices j and k run over the 12 squark mass eigenstates of the MSSM (taken here to be unmixed, but not 
necessarily degenerate), and for SU(3) C , C A — 3 and C q = 4/3 and I q = 1/2. The functions II appearing in eq. (2.7) 
were given in cqs. (5.6)-(5.9) of rcf. [25], and those in eq. (2.8) were given in eqs. (5.5)-(5.10) of rcf. [27]. Note that 
the dependence on all squared masses is explicit in the arguments of these functions. 

Applying eq. (2.6) to this gives the improved equations for the relation between pole and running squared masses: 

M l- lT Q, M Q, = < +^ng>(M^,M?) + /^ 

+I q £ fig-> (Ml ,MlMl)- C,Re[fig> (M^ , M 2 )]^^ (M^. , M|) 



Qj 



-Re [C A U^ (M? ) + I q £ uf' b ) (M? ,M| fc )] ^fig) (M^Mf)}, 

k 9 

h 2 - 



(2.10) 



M?-iT s M- g = m? + / l [c A n^(M?)+7 g ^n( 1 ' fe )(M?,M|p 



cm 



(2,0) 



(M|) 



E nf b) (Af| , m|; + c,/, 2 n? ' c) (Ml , ) + il E E s ? A) (m 2 , m| , m 2 ^ ) 



Re[c A uf-\M 2 ) +I q J2 ^r\Ml Ml)\ ^_ [c A U^(M 2 ) + I q £ I?< 6) (Mf , M?p 



(2.11) 



Given the running masses and coupling as inputs, the pole masses can now be solved for iteratively. Or, given the 
pole masses and the running gauge coupling as inputs, the running masses can be immediately extracted. 

I have checked analytically that the imaginary parts of the pole masses given in eqs. (2.10) and (2.11) correspond 
exactly to the gluino and squark decay widths calculated at next-to- leading order in [46]. This is a good reason to 
prefer the improved version eqs. (2.10) and (2.11) over eqs. (2.7) and (2.8), and more generally eq. (2.6) over eq. (2.5). 
In section IV, I will present a numerical comparison of thes e equations. 



III. Three-loop contributions to the gluino pole the gluino is a color octet, and so is effectively more 
mass 



Radiative corrections to the mass of the gluino in the 
MSSM are particularly large, for two reasons. First, 
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strongly coupled than a quark in the fundamental rep- 
resentation. Second, it couples to 12 squark/quark pairs. 
If one expresses the gluino pole mass in terms of the run- 
ning mass evaluated at itself in a non-decoupling scheme, 
then the squark-mediated corrections are large and grow 
logarithmically with the ratio of the squark to gluino 
masses. One can exploit this by using effective field the- 
ory and renormalization group methods to obtain the 
logarithmically-enhanced parts. In this section, I will 
use this strategy to evaluate the three-loop gluino pole 
mass in the formal limit M? ;§> M|, neglecting terms 



of three- loop order that are suppressed by M-/M?, 
including all terms of order L 3 , L 2 , and L, where 



but 



ln(M?/M? 



(3.1) 



I will also find the coefficient of the L° term, up to a single 
(presently) unknown and plausibly sub-dominant match- 
ing coefficient. This analysis neglects squark mixing and 
non-degeneracy, electroweak effects, and Standard Model 
quark masses, for simplicity. The same method could 
quite easily be extended to include all terms of order 



a n s L n , a n s L n ~\ and a^L n ~ 



(3.2) 



at arbitrary loop order n, but the residual unknown 
three-loop order contributions arc likely to be larger than 



the contributions from n > 4. In the following, I will use 
the following convenient notations for other logarithms: 



ln(M?/Q 2 ), 
ln(M|/Q 2 ), 
ln(m|/m|), 
ln(m?/Q 2 ), 

->2 



ln^/Q 2 ), 



(3.3) 
(3.4) 
(3.5) 

(3-6) 
(3.7) 



where Q is the renormalization group scale, and the run- 
ning gluino and squark masses in the last three definitions 
are taken to be in the full theory (with squarks included) . 

The starting point is the three-loop pole mass for a 
color octet Majorana fermion (g) in the presence of 6 light 
quark flavors (u,d,s,c,b,t). This is the same fermion 
content of SUSYQCD as found in ref. [48] and the ex- 
treme limit of "split supersymmctry" [49]; it is the effec- 
tive theory in which squarks have been decoupled. The 
three-loop gluino pole mass in this model is almost known 
from ref. [50], up to a single ambiguity that is resolved 
in the Appendix of the present paper. The result can be 
written as: 



J 



Mg = m^Q^l+^a^ig + a^] +£ 2 ^ (3.8) 



r 



where M g is the gluino pole mass, and the hats on the 
symbols m g and h = gl/16n 2 are used to distinguish the 
running MS parameters in the effective theory without 
squarks. The coefficients appearing here are: 



,(3) 



= -9, 
36tt 2 In 2 



54£ 



12, 

339 + 33tt 

-282, 4 2) = 63 

172399/18 + 73367tt 2 /10 

+12330C3 - 121tt 4 /2 
+552 In 4 2 + 13248Li 4 (l/2) 

+675^ 2 C 3 - 1890C 5 , 



(3.9) 
(3.10) 
(3.11) 
12672tt 2 In 2 



48tt 2 In 2 2 



(3.12) 



,(3) - 



= -89 32 - 6 2 77r 2 + 6847r 2 ln2 - 306C 3 , (3.13) 



,(3) 



3093, 



4 3) = 

section, 



399. 



(3.14) 



Note that in this section, the gluino pole mass is 
real, because the gluino has no allowed decays in pure 
SUSYQCD in the case that all squarks are heavier. 

The terms in eq. (3.8) that depend on L g , and thus ex- 
plicitly involve the renormalization scale Q, are obtained 



from the renormalization group equations for the gauge 
coupling and running gluino mass in the effective theory: 



dQ 



-h = 



h 2 b^ +h 3 bi 2) +h 4 b^ + . 



dQ 
where 



uit, = m, 



Ml 



4," 

1,(3) 



1.(1) 



1.(3) 



1.(2) 



= -io, 

= 4168/3, 
= -18, b) 
= -2774+ 1440C 3 - 



44, 



J2) 



-228, 



(3.15) 
(3.16) 

(3.17) 
(3.18) 
(3.19) 
(3.20) 



The last coefficient does not seem to be obtainable di- 
rectly from the results in the published literature, which 
only deal with theories with a single type of fermion rep- 
resentation. However, it can be inferred from an unpub- 
lished paper of O.V. Tarasov [51]. This is explained in 
the Appendix of the present paper. 
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The results given above are not what is needed in or- 
der to discern an organizing principle for supersymmctry 
breaking. Instead, one needs to obtain the running pa- 
rameters in the full theory with squarks included. This 
can be achieved from matching conditions between the 
effective theory parameters (h,fh g ) and the full theory 
parameters (h, m§, rr&). In these equations, the effective 

theory parameters (h,rh g ) are always in MS, while the 
full theory parameters (h, trig, m 2 ^) are always in the DR 

scheme here. (It may also be possible to use DR for the 
effective theory, but then there are subtleties associated 
with evanescent couplings in non-supersymmetric theo- 
ries [52].) Now one can take the gluino pole mass from 
cq. (3.8) and use the matching to rewrite it in terms of 
the parameters of the full theory. In the approximation 
used here, non-renormalizable terms are not included in 
the effective field theory Lagrangian, which corresponds 
to neglecting all contributions suppressed by powers of 
M|/M?. Experience shows that expansions of radiative 
corrections in such mass ratios typically converge quite 
quickly, so I suspect that it is reasonable to hope that the 
results below will be approximately valid even if the typ- 
ical squark mass is not very much larger than the gluino 
mass. 

The renormalization group equations for the parame- 
ters of the non-decoupled theory are: 



dQ 



h 



h 2 b^ + h 3 b? + h 3 bf 



+ 



(3.21) 



d_ 

dQ' 
d 



Q—m g 



= m g [hbV+h 2 b%\+hH£\ + ..] (3.22) 



Q dQ m l 



where [53]: 



m 



+m| [h*b™ + .. 



(3.23) 



b£ ] = -6, b {2) = 28, b (3) = 694/3, (3.24) 



and [54, 55]: 



&L n) - nb\ 



h ) 



(3.25) 



and [16, 56]: 



b {1 \ 

,(2h) 



-32/3, 
64. 



k (2o) 



-128/3, 



(3.26) 
(3.27) 



The two-loop gauge coupling matching condition has 
been obtained in ref. [57]: 



h = h 1 + /i(-1 + 2Lq) + h 2 (23- 12Lq + 4L|) 
For the gaugino mass matching condition, I find: 

1 + /i(6Lq) + ^ 2 (14tt 2 /3 - 176 + 133^ + 18L|) 
+h\cZ\ + [1713 - 720C 3 ]£q + 615L| + 48L|) + . . . 



(3.28) 



m 5 = m 5 



(3.29) 



I obtained the two-loop coefficients here by a direct comparison of the two-loop gluino pole squared mass as found in 
the full theory with DR and in the effective theory with MS, using the results of ref. [27]. The three-loop logarithmic 
terms in eq. (3.29) are obtained from the renormalization group equations (3.15)-(3.27). Unfortunately, the coefficient 
Cm s remains unknown, and seems quite difficult to calculate. 

The preceding results can now be used straightforwardly to obtain the three-loop result for the gluino pole mass in 
terms of the non-decoupled theory parameters. I find: 



3] 



M- 9 = m- g (Q) (l + h[a$ + c$L- g + a^L^ +h 2 [a$ + a$>L s + o$Lq + c$L 2 g + a^L- g L^ + a$L 
+n [a 0Q + a 1Q L g + a 01 l, q + a 20 L, g + a n ^ g ^ Q + a Q2 l,^ + a 30 L, g + a 21 ^ g ^ Q + a 12 ^ 9 ^q + « 3 ^q)} 
where 



^00 

-l {2) 

no 

(3) 
2 00 



= 12, 



no 



-9, 



>(!) - 



*01 



-273, 

.(3) 



,(2) 



229, 



= 6, 

,(2) - 



(2) 
x 00 



151 + 113tt 2 /3 - 36tt 2 In 2 + 54C 3 , 



63, 



,(2) 



-72, 



,(2) 
Z 02 



18, 



(3.30) 

(3.31) 
(3.32) 



127147/18 + 73267tt 2 /10 - 12600tt 2 In 2 + 12222< 3 - 121tt 4 /2 - 48tt 2 In 2 2 
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+552 In 4 2 + 13248Li 4 (l/2) + 675tt 2 C 3 - 1890Cs, (3.33) 

a$ = 6483 + 330tt 2 - 3607r 2 ln2 - I8OC3, 4o = -6991 - 669tt 2 + 684tt 2 In 2 - 306C 3 , (3.34) 

4o = 2967, off = -3855, a$ = 1023, 4o = -399, (3.35) 

4? = 630, 4f = -306, 4s = 48 - ( 3 - 36 ) 



It is also possible to rewrite this result so that the logarithms involve running masses, by using the two-loop relation 
between the squark running and pole masses in the formal limit M? » M|, obtained from rcf. [25] using eq. (2.10) 

of the present paper: 

/932 128tt 2 



M l- lT Q M Q = m2 Q + M l 



h\(l in) + h 2 {^ + ^-(1 + 8 ln2) + 24L$ - nr(^ - ^- - &L Q ) }] . (3.37) 

The result, formally equivalent to eq. (3.30) up to terms of four-loop order, is: 

M- g = m- g (Q) (l + h[a$ + c$l- g + o$Iq] +h 2 [a$ + a^i- g + 4l% + c$l% + a^£ g e Q + 

+h 3 [4 3 } + a$£~ g + a^i Q + agk 2 + a^l g £ Q + nfflfy + a <$t§ + + 4fV| + ) . (3-38) 

where the new coefficients are: 

a$ = -49 + 1137r 2 /3 - 367r 2 ln2 + 54C 3 , 4o = - m > 4? = 121 i ( 3 - 39 ) 

a$ = c£l + 60895/18+ 199541tt 2 /30 - (35792/3)7r 2 ln2 + 11250C 3 - 121tt 4 /2 - 48tt 2 In 2 2 

+552 In 4 2 + 13248Li 4 (l/2) + 675tt 2 C 3 - 1890C 5 , (3.40) 

= 809 - 6697r 2 + 6847r 2 ln2 - 306C 3 , 4? = 837 + 330tt 2 - 360tt 2 In 2 - 180<C 3 , (3.41) 

4o = 294, = -723, 0^2 = 159 - (3-42) 

For practical calculations, it is useful to extract from the above expressions the contributions that can be consistently 
added to the complete two-loop results of eqs. (2.8) and (2.11), or more generally eq. (2.5) or (2.6). From eq. (3.38), 
I obtain: 

A cq (2 5) or (2 8)^ = ^"^ + dwi + ^ + ^ + dmt ' a + dllU ' 9 + 

+d 02 £ 2 g + d 12 U 2 g + d 03 e g ] (3.43) 

with coefficients: 

d oa = 2c^l + 50311/9 + 213101tt 2 /15 - (74176/3)tt 2 In 2 + 23796<C 3 - 121tt 4 - 96tt 2 In 2 2 

+ 11041n 4 2 + 26496Li 4 (l/2) + 13507r 2 C 3 - 3780C 5 , (3.44) 

d w = 3990+ 11127r 2 - 11527r 2 ln2 + 288C 3 , d 20 = 2 2 02, d 30 = 3 1 2, (3.45) 

doi = 3826 - 9047r 2 + 8647r 2 ln2 - 1296C 3 , d u = -2280, d 2 i = -864, (3.46) 

d 02 = -384, = 648, d 03 = -108, (3.47) 

and from eq. (3.30): 

A eq (2 6) or (2 11) M § = h * M t [ e °° + &wL + e2oIj2 + + e ° lL ~ a + eilLL ' a + e ° 2 ^] ' ( 3 ' 48 ) 

where 

e 00 = 2c^l + 91507/9 + 59707tt 2 /5 - 22608tt 2 In 2 + 20556Cs + 1104 In 4 2 - 96tt 2 In 2 2 

-121tt 4 + 26496Li 4 (l/2) + 1350tt 2 C 3 - 3780C 5 , (3.49) 

eio = 1410 - 6967r 2 + 5767r 2 ln2 - 2304C 3 , (3.50) 

e 20 = -2310, e 30 = 312, e m = -314, e n = -504, e 02 = 126. (3.51) 

[The coefficients of L 2 L g and LL 2 ~ and L g in eq. (3.48) vanish.] Writing these in a convenient numerical form, I find: 

A (3) ,„_,. , 00 +4 2 = a|m? [0.00101c^ +9.616 + 3.744^ + 1.110^ + 0.157^- 0.375£- q 
eq. (2.5) or (2.8) 9 6 91 m » 9 
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-1.149« § - 0A35£ 2 £g - 0.194^ + 0.327«? - 0.0544^?] , (3.52) 

A (3) , oc , ,oii\ M ? = a|Mf [O.OOlOlc^ + 5.992 - 2.161L - 1.164L 2 + 0.157L 3 
eq. (2.6) or (2.11) s t> g l m g 

+0.158Lg - 0.254LLg + 0.0635L 2 ,] . (3.53) 

When applied to realistic models, the two-loop gluino pole mass should be calculated including all relevant effects 
including squark mixing and Yukawa and electroweak couplings using eq. (2.5) or (2.6), and then the appropriate 
corresponding formula eq. (3.43) or (3.48) can be added with an approximate overall squark mass scale parameterized 
by either £ or L. (However, the three-loop contributions found here must be eschewed if most of the squarks are not 
heavier than the gluino.) 

The parts of eqs. (3.52)-(3.53) that do not contain logarithms came from two sources: the non-logarithmic part of 
eq. (3.8), and the unknown three-loop gluino mass matching coefficient c m ~ defined in eq. (3.29). It is useful to note 
that the dominant part of the contribution from eq. (3.8) is due to loop diagrams containing only gluons and light 
quark internal lines. Furthermore, there is a significant partial cancellation in this non-logarithmic piece, not due 
to supersymmetry, which is not present in the effective theory, but to the accident of the number of quarks in the 
Standard Model. To see this, one can tag the contributions according to the number n of closed gluino loops in each 
diagram, and also treat the number of light quarks uq (equal to 6 in the real world) as a variable, using eq. (A. 15). 
Then, in eqs. (3.52) and (3.53) respectively: 

9.616 -> 33.312 - 26.634(n Q /6) + 3.489(n Q /6) 2 - 0.027(n Q /6) 3 + [-0.598 - 0.0095(n Q /6)]n + 0.084n 2 , (3.54) 
5.992 -> 16.703 - 10.543(n Q /6) + 0.667(n Q /6) 2 + 0.027(n Q /6) 3 + [-1.072 + 0.126(n Q /6)]n + 0.084n 2 . (3.55) 

Thus the diagrams with no closed heavy particle loops dominate the final result, and would even more so if it were not 
for the tendency of gluon and quark loops to cancel. In the numerical studies of the following section, I will simply 
neglect the effects of the unknown coefficient c m '~ , since it comes from diagrams with at least one heavy squark loop 

and therefore is plausibly less significant than the other non-logarithmic contributions. Note that \dm s \ would have 
to exceed 10 4 in order for it to contribute 1% to the gluino pole mass formula. 

It is also useful to observe that there will typically be a significant cancellation between the logarithmic and non- 
logarithmic contributions in eq. (3.53). Therefore, that version of the three-loop contribution to the pole mass is 
actually considerably smaller than one might have naively suspected. 

IV. Numerical results 

For purposes of illustration, consider a simplified model, with all squarks degenerate and unmixed, and quark 
masses and electroweak effects neglected, as in the previous section. The pertinent Lagrangian parameters are then 
the running SUSYQCD coupling as(Q) and the gluino and common squark masses in the DR scheme. Throughout 
this section, I will fix as{Mg) = 0.095. 

In the left panel of Figure 1, I compare different computations of the ratio of the real part of the gluino pole mass 
to the running mass evaluated at the pole mass, Mg/m,g(Mg). The two-loop computations of eqs. (2.8) and (2.11) for 
the pole mass agree to better than 1% for Mq/M 9 < 1.55, but the disagreement increases for larger values of that 
ratio, and reaches 2.7% when Mq/M 3 — 3.5. It is in just this regime that the three-loop contributions to the gluino 
pole mass found above should be reliable and important; I will return to this below. 

In the right panel of Figure 1, the same comparison is done for the ratio of the real part of the squark pole mass to 
the running mass. Here, the different two-loop computations are in extremely close agreement over the entire range. 
Furthermore, the overall magnitude of the radiative corrections is much smaller than for the gluino. I conclude that 
purely theoretical uncertainties for squark masses are probably under control at a level of much better than 1 per 
cent. (The steep "cliff" at the left side of the graph reflects the fact that a much heavier gluino makes a large negative 
radiative contribution to the squark pole squared mass.) 

I expect that the solid lines in Figure 1, reflecting the calculations of eqs. (2.10) and (2.11), are more reliable 
than those of eqs. (2.7) and (2.8). As discussed above, one reason for this expectation is the fact that the former 
equations do a much better job of approximating the decay widths (the imaginary parts of the pole mass) in the 
near-threshold region. This is illustrated in Figure 2. First, the gluino width T g as calculated from the pole mass 
using eq. (2.8) actually vanishes for all Mq > 0.969M§, rather than for Mq > Mg as dictated by kinematics. The 
reason for this failure is that the width in the pole mass derives from the imaginary parts of loop integrals which, in 
that approximation, depend on running masses in the propagators instead of physical masses. The approximation of 
eq. (2.11) does not have this problem, and exactly reproduces the direct next-to- leading order width calculation of 
ref. [46]. 
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FIG. 1: The ratios of the gluino and squark pole masses to the running masses, Mg/rrig(Mg) in the left panel and Mq/itIq^Mq) 
in the right panel, as functions of the ratio of the common squark pole mass to the gluino pole mass, Mq /M§ . The different lines 
correspond to the two-loop results of eqs. (2.8) and (2.8) and the one-loop truncations of the same formulas. For simplicity, 
here all squarks are taken to be degenerate and unmixed, and quark masses and electroweak effects are neglected. The 
renormalization scale use for the computation is Q — Mg for the left panel and Q = Mq in the right panel, and the SUSYQCD 
coupling is fixed to as(Mg) — 0.095 in both cases. 




FIG. 2: The gluino and squark widths as extracted from the complex pole masses, Tg/Mg in the left panel and Vq/Mq in the 
right panel, as functions of the ratio of the squark pole mass to the gluino pole mass, in the near-threshold region of parameter 
space. The different lines correspond to different two-loop approximations, as defined in the text. The widths as computed 
using eqs. (2.10), (2.11) agree exactly with the direct next-to-leading order width calculation of ref. [46], while the widths as 
computed using eqs. (2.7), (2.8) fail to agree with the decay kinematics dictated by the physical masses near threshold. 



A similar situation holds for the squark width, as shown in the right panel of Figure 2. In fact, the calculation of 
eq. (2.7) gives slightly negative (and therefore unphysical) values for the width in a narrow range on either side of 
the physical threshold. This is because the two-loop pole mass contribution to the width overcompensates for the 
one-loop contribution, which is strictly positive for all Mq > 0.969M§. Again, the two-loop calculation of eq. (2.10) 
gets the kinematics correct, and precisely reproduces the next-to-leading order calculation of ref. [46]. 

I next turn to the effect of the partial three-loop contributions, derived in section III, on the gluino pole mass. 
This is shown in Figure 3. Strictly speaking, the three-loop calculations given here are only valid in the formal limit 
M? ^> M? but in any case the applied correction is small for squark masses just above the gluino mass, so I have 
taken the liberty of showing the entire range Ma > M g . The two three-loop approximations are much closer to each 



10 




FIG. 3: The ratio of the gluino pole mass to the running mass, Mg/mg(Mg), as a function of the ratio of the squark pole masses 
to the gluino pole mass, Mg/Mg, as in Figure 1. The different lines correspond to different two- loop and partial three- loop 
approximations, as defined in the text. 



2-loop, eq. 
2-loop, eq. 



(2.8) 
(2.11) 



1 



Q/M 



1.28 



1.27 



1.26 



,1.25 



1.24 



1.23 



1.22 




■ 2-loop, eq. (2.8) 
■■ 2-loop, eq. (2.11) 

■ 3-loop, eqs. (2.8), (3.52) 
- 3-loop, eqs. (2.11), (3.53) 



1 



gluino 



Q/M 



gluino 



FIG. 4: Renormalization-group scale (Q) dependence of the ratio of the gluino pole mass to the running mass, Mg/rrig(Mg). The 
left panel shows the two-loop approximations of eqs. (2.8) and (2.11), for Mg/Mg = 0.9 (lower pair of lines) and Mg/Mg = 1.5 
(upper pair of lines). The right panel shows the same information, and also the three-loop approximate contributions of 
eqs. (3.52) and (3.53), for M- a /Mg = 3. 



other than the corresponding two-loop approximations, just as one might have hoped. They differ by less than 1% 
even for Mg/Mg — 5. It is also noteworthy that both three-loop results are closer to the two- loop approximation of 
eq. (2.11) than they are to eq. (2.8), providing some circumstantial evidence for the superiority of eq. (2.11). 

Finally, consider the renormalization group scale dependence of the calculated relationship between the pole mass 
and the running mass of the gluino. Numerical results are shown in Figure 4, for three ratios Mg/Mg = 0.9, 1.5, and 
3. In each case, the ratio of the pole mass Mg to the running mass evaluated at the pole mass, rrig(Mg) is computed 
for a fixed model in terms of the renormalization scale Q at which the calculation of the pole mass is performed. The 
renormalization group equations (3.21)-(3.27) are used to run the running parameters between different values of Q. 
Comparing the two-loop results, the approximation of eq. (2.11) is slightly more stable than that found using eq. (2.8), 
although both are quite acceptably scale-invariant for Mg/Mg = 0.9 and 1.5. In the case of M g /Mg = 3, shown in the 
right panel, I also include the three-loop contributions of eqs. (3.52) and (3.53). They exhibit a still further improved 
scale dependence; this is encouraging but cannot be counted as a surprising triumph, since the explicit Q dependence 
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of the three-loop contribution to the pole mass came from n othing other than the three-loop beta functions. 



V. Conclusion 

In this paper, I have argued in favor of a reformulation 
of the two-loop approximation between pole and running 
squared masses. As an improvement over eq. (2.5), equa- 
tion (2.6) has general applicability. It was applied here to 
the specific case of gluinos and squarks in the SUSYQCD 
sector of the MSSM in eqs. (2.10) and (2.11). I also used 
the method of effective field theory to obtain a partial 
three-loop approximation to the gluino pole mass, when 
squarks are heavier. The agreement between the three- 
loop gluino pole mass results and the two-loop approxi- 
mation of eq. (2.11) provides evidence that the method 
of expanding running masses about the real part of pole 
masses in the loop corrections provides better accuracy. 
Another piece of evidence in favor of this conjecture is the 
agreement of the imaginary part of the pole mass with a 
direct calculation of the width. Also, the improved renor- 
malization scale dependence is at least consistent with it. 
If the LHC discovers strongly interacting supcrpartners, 
then the quest to decipher the organizing principle be- 
hind supersymmetry breaking should eventually benefit 
from the improved results presented here, as well as sim- 
ilar applications to the rest of the sparticle spectrum. 

Appendix: Three-loop results for gauge theories 
with fermions in arbitrary representations 

In this Appendix, I compile the following results for 
a general gauge theory renormalizcd in the MS scheme 
with fermions in arbitrary representations and no scalar 
fields: 

• the three-loop beta functions for gauge couplings, 

• the three-loop beta functions for fermion masses, 
and 

• the three-loop relation between the pole and run- 
ning MS masses, in the limit that there is only one 
non- vanishing fermion mass parameter. 

Each of these results has appeared before in the case 
of theories that are QCD-like (containing a single gauge 
group and a single type of fermion representation) . How- 
ever, there arc non-trivial ambiguities in inferring the 
three-loop results for general theories from the published 



literature. The purpose here is to resolve these ambigui- 
ties for use in the main text of the present paper and for 
future reference. 

To set notation, consider a theory with a gauge group 
G which is the product of one or more simple or U(l) 
gauge groups G a , each with a distinct MS gauge coupling 
g a . Results will be written in terms of the combinations 

h a ee gl/Wn 2 . (A.l) 
Suppose further that the two-component Weyl fermions 
of the theory transform in possibly distinct representa- 
tions of the gauge group labelled by R. Each Dirac (Ma- 
jorana) fermion consists of two (one) such Weyl fermions, 
so this entails no loss of generality. The quadratic 
Casimir invariants of the adjoint and fermion represen- 
tations of each gauge group are written as C a (A) and 
C a (R), respectively. The normalization is such that 
C a (A) = N for SU(N), and C a (R) = (N 2 - l)/2N 
when R is a fundamental representation of SU (N) . The 
Dynkin index of each representation R is written as 
I a (R), in a normalization such that the fundamental rep- 
resentation of SU(N) has index 1/2 for a Weyl fermion. 
I also define the invariants: 

S a = ^la(R), (A.2) 
R 

Sab = 2/„(i2)C 6 (ii), (A.3) 

R 

Sabc = ^4(ii)C t (iJ)C c (iJ). (A.4) 

R 

Note that a Dirac fermion contributes twice to each of 
these sums. For example, a Dirac fermion in a funda- 
mental representation of SU(N) will contribute 1 to S a , 
and a Dirac fermion with charge q under a U(l) gauge 
group will contribute 2q 2 to the corresponding S a . 

The three-loop beta function for each of the MS gauge 
couplings is: 

Ph - = Q ^Q K = ^ + ^ + ^ + ■ ■ ■ (A ' 5) 
where the terms in the loop expansion are: 

0™ = h 2 a [-22C a (A)+AS a }/3, (A.6) 



Ht] = h 3 a C a (A)[-68C a (A) + 20S a ]/3 + 4h 2 a h b S ab , (A.7) 
Ph! = h i a C a (A){-2857[Ca(A)} 2 + 1415Ca(A)Sa-79[S a } 2 }/27 

+h 3 a h h 8C a (A)Sab + h 2 a h 2 b S ab [l33C b (A) - 22S b ]/9 - 2h 2 a h b h c S abc , (A.8) 
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with the indices b and c implicitly summed over in each 
term where they appear. The special case of this result 
for a QCD-like theory with a single gauge group com- 
ponent and a single type of fermion representation was 
given in ref . [58] . (The four- loop result has also been ob- 
tained in QCD-like theories, in [59]). In that special case, 
the three-loop terms proportional to h^hbC a (A)S a b and 
h 2 a h 2 S a bCb{A) are combined, and the term proportional 
to h 2 h 2 S a bSb could in principle combine with a term pro- 
portional to h^hbSabSa, which is actually absent. These 
two ambiguities have been resolved by considering the 
special case of the electromagnetic coupling beta func- 
tion with QCD effects included; see for example eq. (54) 
of [60] and eqs. (9)-(10) of [61]. _ 

Next consider the three-loop MS beta function for the 
mass m of a fermion transforming in a representation F, 



given for a QCD-like theory in [51]. Chiral symmetry 
guarantees that if there are several such masses, they 
each run independently. As far as I know, the result for 
a theory with different fermion representations has not 
been given directly in the published literature, but can 
be inferred by considering the results given for individual 
classes of diagrams in [51]. The result is: 



dQ 



where the terms in the loop expansion are: 
tf» = -6h a C a (F), 



(A.9) 



(A.10) 



= h 2 a C a (F)[-97C a (G) + 10SJ/3 - 3h a h b C a (F)C b (F) 

/1 3) = hlC a {F){- l -^[C a {A)] 2 
129 



(A.11) 



556 
~27 



48C 3 



C a (A)S a + 7 ^(S a f} 



+h 2 a h b C a (F){ — C a (A)C b (F) + S a Cb(F) + [45 - 48( 3 }S ab } - 129h a h b h c C a (F)C b (F)C c (F). (A.12) 



with indices a, 6, c summed over in terms in which they 
appear. The terms proportional to C a (F)S a Cb{F) and 
C a (F)S a b are combined in the case of quark masses in 
QCD, and it is this ambiguity that has been removed 
using the results inferred from [51]. (The QCD-like case 
has been extended to four-loop order in [62].) 

Finally, consider the three-loop fermion pole mass. Let 
the two-component fermions consist of massless fermion 
species with representations labelled by r, as well as de- 
generate massive fermion(s) with representation labelled 
by F and a running mass m{Q). This is only technically 
natural if F is irreducible (as for Majorana fermions), or 
consists of an irreducible representation and its conjugate 
(as for Dirac fermions), or if F consists of three or more 
degenerate copies of a single irreducible representation 
and/or its conjugate (a situation for which I know of no 
examples in proposed extensions of the Standard Model). 
Therefore, it is assumed here that all of the irreducible 
representations labelled by F have the same Casimir in- 
variant C a {F) and index I a {F). The invariants S a , S a b 
previously defined are now separated into contributions 



I 

from the massless and massive fermions: 
S L a = 5>(r), S« = £/ (F), (A.13) 

r F 

SL = £>(r)C&(r), S5 = 2j„(F)C 6 (F). (A.14) 

r F 

(Again one must remember that the representations 
are defined for two-component fermions, so each Dirac 
fermion contributes twice to the appropriate sums.) 
Then the fermion pole mass M is related to the running 
MS mass m evaluated at a renormalization scale Q = M 
by: 

m(M) = M \l + x (1) + x (2) + x (3) +...], (A.15) 

where the loop expansion terms are: 

x {1) = -h a AC a {F), (A.16) 



,(2) 



,0) 



h 2 a Ca(F){C a (A)[-llll/24 + 4ir 2 /3 - 4tt 2 In 2 + 6( 3 ] + S% [71 + 8tt 2 ]/3 + Sf [143 - 16tt 2 ]/3} 
+h a hbCa{F)Cb{F)[%Ti 2 \n2-bn 2 - 12C3 + 7/8], (A.17) 
64^ a C a (^) \h h h c Cb{F)C c {F)df } + h a h b {C a (A)Cb(F)d^ + 2S^ b df A + 2S^C b (F)d ( ^ + 2Sf 6 4 3) } 



+h 2 a {[C a (A)] 2 d^ + 2C a {A)S^ + 2C a (A)Sf $ 
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4[Sf] a 4 8) +4[tf] a d8>} , (A.18) 
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with indices a, 6, c summed over wherever they appear. 
The coefficients d„ were found in ref. [50] for n = 
1,2,3, 5, 6, 7, 8, 9, 10 and will not be repeated here. The 
remaining coefficients and d^g were combined into 

(3) 

a single coefficient d\ in that paper, since those terms 
are indistinguishable in the special case of a single type 
of fermion representation. In this paper, I need the gen- 
eralized result: 



,(3) 
AA 

.(3) 
4B 



827/384 + 71-716 

•2 



tt 4 /240- 11C 3 , 



85 95 o 11 2l „ 11 11 

1 ir 2 ir 2 In 2 H C3 7T 

1152 144 9 3 s3 216 

+Vln 2 2 + -ln 4 2- 



(A.19) 

4 



;LU(l/2). 



(A.20) 



provided by Kirill Melnikov. 

In the application of the present paper, the effective 
theory with squarks decoupled consists of an SU(3) C 
gauge theory with 6 flavors of "massless" Dirac fermion 
quarks and 1 massive color octet Majorana gluino. 
Therefore, h a — as/47r and the relevant group theory 
invariants are: 



C a (A) = 


C a (F) = 3, 






(A.21) 


Sa = 6, 




oL 

^aaa 


= 32/3, 


(A.22) 


Sf = 3, 


S H - 9 

J aa a ; 


oH 

^aaa 


= 27, 


(A.23) 


S a =9, 


Saa = 17, 


Saaa 


= 113/3. 


(A.24) 



I obtained d 



(3) 



AA by a direct computation of the corre- 

(3) 

sponding three-loop diagrams, and then obtained d\g = 

<M — ^4/1 using the result for d\ provided in ref. [50] . 
This was also checked independently using a slight mod- 
ification of the computer code used in ref. [50], kindly 
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